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The spectrum of low-lying eigenvalues of overlap Dirac operator in quenched SU (2) lattice gauge 
theory with tadpole-improved Symanzik action is studied at finite temperatures in the vicinity of 
the confmement-deconnnement phase transition defined by the expectation value of the Polyakov 
line. The value of the chiral condensate obtained from the Banks-Casher relation is found to drop 
down rapidly at T = T c , though not going to zero. At T' c w 1.5 T c ~ 480 MeV the chiral condensate 
decreases rapidly one again and becomes either very small or zero. At T < T c the distributions 
of small eigenvalues are universal and are well described by chiral orthogonal ensemble of random 
matrices. In the temperature range above T c where both the chiral condensate and the expectation 
value of the Polyakov line are nonzero the distributions of small eigenvalues are not universal. 
Here the eigenvalue spectrum is better described by a phenomenological model of dilute instanton 
- anti-instanton gas. 



PACS numbers: 12.38.Gc; 25.75.Nq 

I. INTRODUCTION 

It is well known that in the limit of zero quark masses 
the classical action of QCD is invariant under the ex- 
change of quark fields with different chirality. This clas- 
sical chiral symmetry, however, appears to be broken in 
quantum theory. Its spontaneous breakdown gives rise 
to massless Goldstone particles (pions), and to nonzero 
expectation value of the chiral condensate ( rptp ) . Since 
pions are the only massless hadrons in this limit, one 
can describe infrared QCD using an effective chiral the- 
ory where pions are the only degrees of freedom [1, 2]. 
At sufficiently high temperatures chiral symmetry is re- 
stored again, and the effective chiral theory is not valid 
any more, which means that at such temperatures the ef- 
fective degrees of freedom arc not pions. It is commonly 
believed that the restoration of chiral symmetry is as- 
sociated with the confinement-deconfinement transition, 
with the chiral condensate being the order parameter. 
For quenched gauge theories the chiral condensate is not 
an exact order parameter, since in this case it does not 
correspond to expectation value of any physical quantity 
and can be defined only indirectly, via the Banks-Casher 
relation [31: 



SE|(#)| = lim lim ^ 

A— >0 V— >oo V 



(1) 



where p (A) = S (A — Aj) ) is the density of eigenvalues 

i 

Aj of the Dirac operator and V is the volume of the four- 
dimensional box in which the theory is considered. Thus 
the equation (1) relates the chiral condensate and the 
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density of small eigenvalues of the Dirac operator. 

In quenched SU (N) gauge theories a commonly used 
order parameter is the expectation value of the Polyakov 
line, which is associated with Zjq center symmetry and 
which is equal to zero in the confinement phase, where 
center symmetry is unbroken [4] . Thus the Polyakov loop 
and the chiral condensate are the proper order parame- 
ters for the deconfinement phase transition in the lim- 
its of infinite and zero quark masses, respectively. How- 
ever, numerical simulations indicate that both the chiral 
condensate and the Polyakov loop can be used as ap- 
proximate order parameters even when the correspond- 
ing symmetries are absent. In particular, for QCD with 
dynamical quarks Polyakov line changes rapidly in the 
vicinity of the deconfinement phase transition [5]. The 
measurements of the chiral condensate using the stag- 
gered Dirac operator in quenched SU (2) gauge theory 
has also shown that it goes to zero in the deconfine- 
ment phase [(>]. However, in some recent lattice studies of 
quenched SU (2) lattice gauge theory with chirally invari- 
ant lattice Dirac operator it was found that even above 
the deconfinement phase transition there arc still some 
small eigenvalues Xi which yield nonzero chiral conden- 
sate in the Banks-Casher relation (1) [7, 8]. To be more 
precise, it turns out that in the deconfinement phase with 
spontaneously broken center symmetry chiral condensate 
behaves differently depending on the sign of the Polyakov 
loop. For field configurations with positive Polyakov loop 
the chiral condensate eventually goes to zero at some 
temperature above T c , while for configurations with neg- 
ative Polyakov loops the chiral condensate stays nonzero 
at all temperatures which were considered [8]. Since for 
SU (2) gauge theory with dynamical fcrmions at temper- 
atures above the deconfinement phase transition positive 
values of Polyakov loop are favoured, at T > T c it seems 
reasonable to consider only field configurations with pos- 
itive Polyakov loops in order to capture this important 
feature of the full QCD already in the quenched approx- 
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imation [7, 8]. 

The value of the chiral condensate is usually extracted 
from the results of lattice simulations by measuring the 
probability distribution p (A) of the eigenvalues of the 
Dirac operator V = 7^ (<9 M — iA^) and by using the 
Banks-Casher relation (1). A much more precise method 
was proposed in [9, 10]. This method is based on the fact 
that the properly rescaled eigenvalue distribution p (A) 
at A — > is universal and coincides with the eigenvalue 
distribution of the so-called chiral ensemble of random 
matrices [2, 11]: 



y 1 ™ EV K^) = Ps {z) 



(2) 



The form of the function ps (z) depends only on the 
global symmetries of the Dirac operator and on the topo- 
logical charge of the gauge fields. One can also con- 
sider the probability distribution p (X m in) of the lowest 
nonzero eigenvalue of the Dirac operator, which is also 
universal if X m i n is rescaled by the factor EV. Numeri- 
cally it was found that for the low-lying eigenvalues the 
universality holds with a good precision already at not 
very large lattice volumes, which can be easily achieved 
in lattice simulations [9], In this case for p(X) with A 
of order of (EV) 1 one has the following approximate 
equalities: 



p{\) = ^Ps (EVA) 

p(X m in) = PS (EVA m i„) 



(3) 
(4) 



Thus in order to find the value of E one can fit the low- 
est edge of the numerically obtained distributions p (A) 
or p(X m in) with the universal functions Aps{cX) or 
Aps (cX m in), where A and c are fitting parameters. The 
parameter c is related to E as E = cV . Since such nu- 
merical procedure involves fitting lattice data with some 
function rather than finding a single value p (0), the pre- 
cision of this method is usually significantly higher than 
that of the Banks-Casher relation. Finite- volume correc- 
tions are also believed to be smaller, which allows one to 
measure the chiral condensate using such small lattices 
as 4 4 [9, 10]. 

In this paper we report on our measurements of the 
spectrum of low-lying eigenvalues of lattice Dirac op- 
erator in quenched SU (2) lattice gauge theory with 
tadpole-improved Symanzik action (see, e.g., the ex- 
pression (1) in [8]) at finite temperatures in the vicin- 
ity of the confincment-deconfincment phase transition 
defined by the expectation value of the Polyakov line. 
In Section II we study the temperature dependence of 
the chiral condensate using the Banks-Casher relation 
[■'!]. We demonstrate that the chiral condensate changes 
rapidly in the vicinity of the deconfincmcnt phase tran- 
sition, though not going to zero. Another interesting 
effect is a second rapid decrease of the condensate at 
T w 1.5 T c w 480 MeV. In Section III we check the 
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TABLE I: Lattice parameters which were used for simulations. 



applicability of the chiral Random Matrix Theory [2] to 
the spectrum of low-lying Dirac eigenvalues at finite tem- 
peratures. It is found that at T < T c the spectrum is 
universal and the Random Matrix theory can be applied, 
while above T c some other model should be used. 

Overlap Dirac operator for quarks in the fundamen- 
tal representation is used to find the eigenvalue density 
p(X) in (3). Advantages of the overlap Dirac operator 
are explicit chiral symmetry at all lattice spacings and 
the existence of exact zero modes for field configurations 
with nonzero topological charge [12]. The parameters of 
our lattice simulations are summarized in Table I. We 
mostly use the same lattice configurations as in [8], thus 
our work can be considered as an extension of the anal- 
ysis performed in that work. Specified numbers of lat- 
tice configurations are the total numbers of field config- 
urations with all topological charges. Following [7, 8], 
at temperatures above T c we consider only gauge field 
configurations with positive expectation values of the 
Polyakov line. For quenched SU (2) gauge theory with 
tadpole-improved Symanzik action recent measurements 
give (3 C = 3.248(2) for Nt = 6, which corresponds to 
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T c = 313.(3) MeV [13]. 



II. TEMPERATURE DEPENDENCE OF THE 
CHIRAL CONDENSATE 

In this Section we present technical details of our mea- 
surements of the chiral condensate and discuss its tem- 
perature dependence as well as finite-volume and finite- 
spacing effects. 

Our results as well as the results of the previous works 
on the spectrum of overlap Dirac operator [7, 14] sug- 
gest that in the vicinity of T c m 320 MeV the spec- 
trum starts developing a sort of "plateau" for A in the 
range 100. . AOOMeV, which becomes wider and lower 
at higher temperatures. At the same time a character- 
istic peak gradually emerges near the origin, below ap- 
proximately 100 MeV. A similar feature of the spectrum 
of small eigenvalues in the deconfinemcnt phase was ob- 
served in [7] for quenched SU (2) and SU (3) gauge the- 
ories with Wilson action and with both staggered and 
overlap Dirac operators and in [8] for quenched SU (2) 
gauge theory with tadpole-improved Symanzik action 
and with overlap Dirac operator. The spectra of 30 and 
50 lowest Dirac eigenvalues for 16 3 x 6 lattice at /3 = 3.40, 
T = 401 MeV and for 20 3 x 6 lattice at /? = 3.325, 
T = 346 MeV, correspondingly, are plotted on Fig. 1 
as an example. 

For our data the height of the peak (as measured from 
the origin) remains practically constant at different lat- 
tice volumes and lattice spacings. Moreover, a close in- 
spection of our results as well as of the results presented 
in [7, 14] suggests that the width of this peak goes to zero 
in the continuum limit a — > 0. In order to demonstrate 
this, we have fitted the spectra of the Dirac operator at 
A < 200 MeV by the function C + Acxp (-BX), where 
the parameter B was considered as the inverse of the 
peak width and C was the value of spectral density at 
the "plateau" (see Fig. 1). The width of the peak ex- 
tracted in this way is plotted on Fig. 3 as a function of 
lattice spacing for 16 3 x 6, 20 3 x 6 and 24 3 x 6 lattices. 
It can be seen that the peak width rapidly drops down 
when the spacing is sufficiently small and the correspond- 
ing temperature is in the vicinity of the deconfinemcnt 
phase transition, and approaches zero as lattice spacing 
becomes smaller. The widths of the peaks at different lat- 
tice sizes agree within error range. We have also analyzed 
the same data as the functions of four-dimensional and 
three-dimensional lattice volumes in physical units, but 
in these cases the agreement between different lattices 
was much worse. Thus it seems that this peak is indeed 
a lattice artifact, since its width goes to zero while its 
height remains practically constant, and in the contin- 
uum limit Dirac eigenmodes which constitute the peak 
turn into exact zero modes. 

Although it has been suggested in [14] that the height 
of the peak diverges, our results in fact do not con- 
tradict those reported in [14]. In [14] the shape of 



the peak was approximated by a function of the form 
p{E) = a(V) E- b , where E = Xa is the virtuality 
A in lattice units, a (V) is some coefficient which in- 
creases with lattice volume V and the optimal value of b 
is b « 0.80. The dependence of the peak parameters on 
lattice spacing was not investigated. If one rewrites the 
expression for p (E) obtained in [14] in physical units, one 
obtains p (A) ~ a (V) a 1_h A~ b . Thus for any finite value 
of lattice spacing the spectral density indeed diverges at 
A — ► 0. However, since b < 1, this divergence and the 
peak near zero disappear in the continuum limit a — > 0. 
The model of dilute instanton - anti-instanton gas used in 
I] to describe this peak also does not involve any phys- 
ical scale, which again indicates that this feature of the 
Dirac spectrum might be a lattice artifact. 

Taking into account the above considerations, we have 
assumed that at finite lattice spacing the spectral den- 
sity at A — > in the Banks-Casher relation (1) should 
be replaced by the best-fitting value, C, of the rescalcd 
spectral density irp (A) jV at the "plateau" . These values 
of C are shown on Fig. 1 with solid horizontal lines. 

At temperatures below 400 MeV we have calculated 
C by averaging the spectral density over A in the 
range 100 . . . 200 MeV, which is a typical extent of the 
"plateau" for such temperatures (see Fig. 1, right plot). 
For higher temperatures the "plateau" becomes signifi- 
cantly lower and wider. Due to limited statistics, this 
"plateau" on the histograms of p (A) typically consists 
of only a small number of randomly distributed individ- 
ual eigenvalues (see Fig. 1, left plot), and the number 
of eigenvalues between 100 MeV and 200 MeV strongly 
fluctuates. In order to reduce these fluctuations, for 
higher temperatures we have extended the averaging in- 
terval up to Amax = 400 M eV. After such an extension 
the values of the chiral condensate above T = 400 MeV 
agree very well for different lattice volumes. This is an 
additional argument in favor of our method to measure 
the chiral condensate at finite lattice spacing. 

The resulting temperature dependence of £ is shown 
on Fig. 2. A comparison of the data for different lat- 
tice volumes shows that the finite- volume effects are sig- 
nificant for 8 3 x 4 and 10 3 x 4 lattices. Starting from 
12 3 x 6 lattice, for T < T c « 320 MeV the chiral con- 
densate stabilizes at S 1 / 3 w 300 MeV. At T w T c the 
condensate rapidly changes from S 1 / 3 « 300 MeV to 
E 1 / 3 w 120 MeV and then practically does not change 
in the temperature range T c < T < 1.5T C ps 480 MeV. 
This value of £ is also practically volume-independent 
for 12 3 x 6, 16 3 x 6, 20 3 x 6 and 24 3 x 6 lattices, which 
suggests that nonzero chiral condensate above T c is not 
a finite-volume effect. Finally, at T w 1.5 T c w 480 MeV 
£ quickly decreases once again. Although the chiral con- 
densate is equal to zero for the highest temperature that 
we have considered (T = 494 MeV, the last point on the 
right on Fig. 2), this can be just the effect of limited 
statistics, as discussed above. This is an interesting as- 
pect of the theory which deserves further investigations, 
since in quenched theory the chiral symmetry is not a 
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FIG. 1: The rescaled density (np (A) /V) 1/a of eigenvalues of overlap Dirac operator for 16 3 x 6 lattice at f3 = 3.40, T = 401 MeV 
(on the right) and for 20 3 x 6 lattice at (3 = 3.325, T = 346 MeV (on the left). Only the first 30 and 50 lowest eigenvalues were 
taken into account for the right and the left plots, respectively (see Table I). 
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FIG. 2: The values of the chiral condensate E 1 / 3 obtained 
from the Banks-Casher relation (1) at different temperatures. 
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symmetry of the action and strictly speaking there are 
no theoretical reasons to expect that the chiral conden- 
sate should eventually go to zero at some temperature. 
Nevertheless, our results suggest that near T w 480 MeV 
there is at least a sort of crossover transition associated 
with rapid change in the value of the chiral condensate, 
and that the commonly accepted deconfincment picture 
with zero chiral condensate at high temperatures is valid 
at least approximately. 



III. COMPARISON WITH CHIRAL RANDOM 

MATRIX THEORY AND THE MODEL OF 
DILUTE INSTANTON - ANTI-INSTANTON GAS 

In this Section we test the universality of the spectra of 
low-lying Dirac eigenvalues at different temperatures by 
comparing them with the predictions of chiral Random 



Matrix Theory [2]. Non- universal spectra are compared 
with the model of dilute instanton - anti-instanton gas 
used in [7]. 

The chiral ensemble of random matrices which deter- 
mines the universal distributions ps (z) and ps {X m in) for 
overlap Dirac operator with SU (2) gauge group is the 
chiral orthogonal ensemble, an ensemble of real random 
matrices with chiral structure [2, 15]. For such ensemble 
the function ps (z) can be expressed in terms of some 
rather complicated integrals which involve Bcssel func- 
tions [15]. These integrals are rather difficult even for 
numerical integration, thus to obtain the function ps (z) 
we have simply generated 5 • 10 5 chiral random matrices 
of size 50 x 50 for topological charges Q = 0,1,2 and 
used the interpolated distributions of their eigenvalues. 
Further increasing the size of the matrices or the num- 
ber of matrices in the ensemble does not change the result 
within the accuracy of several tenth of percent. The func- 
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TABLE II: The values of the chiral condensate E 1 / 3 extracted 
by fitting the Dirac spectra obtained in lattice simulations 
with the universal eigenvalue distributions of chiral orthogo- 
nal ensembles with different topological charges. The values 
of E 1,/3 obtained from Banks-Casher relation are also shown 
for comparison. The symbols "L" or "A" mean that the distri- 
butions of the lowest eigenvalues p (X m in) or of all eigenvalues 
p(A) were fitted. 



tions ps (z) are known in analytical form for topological 
charges [16]: 

PS {z) = (2 + z) /4 exp (-z/2 - z 2 /8) , Q = 

Ps (z) = z/4exp(-z 2 /8), Q=l (5) 

For Q = 2 the function p$ (z) was found numerically, 
similarly to ps (z). 

As an example, the histograms of the distributions 
p (A) and p (A mi „) at temperatures T = 264 M eV (10 3 x 4 
lattice, T/T c = 0.84) and T = 294 MeV (12 3 x 6 lattice, 
T/T c = 0.94) are plotted on Fig. 4 together with the 
functions ps (EVA) and ^yPs {EV\ m i n ) at the val- 
ues of £ which best fit the lattice data. These values are 
summarized in Table II, where the results obtained using 
the Banks-Casher relation are also included for compar- 
ison. It can be seen that for all plots in Fig. 4 the lat- 
tice data agrees with the universal distributions within 
the range of statistical errors. For 10 3 x 4 and 12 3 x 6 
lattices the fits for different topological sectors give dif- 
ferent results for E, and the Banks-Casher relation also 
yields somewhat higher value of the chiral condensate. 
For 16 3 x 6 lattice the agreement between all measure- 
ments is much better. This indicates that for small lat- 
tices finite- volume effects are rather significant. These 
effects, however, influence only the value of £ but not 
the shape of the distribution of small eigenvalues, which 



remains universal. 

The distributions of eigenvalues for 16 3 x 6 lattice at 
T = 316 MeV, which is just above T c , are plotted on Fig. 
5. It can be seen that these distribution are indeed not 
similar to those predicted by random matrix theory, and 
can be characterized by a strong excess of eigenvalues be- 
low approximately 100 MeV. This characteristic peak in 
the density of the eigenvalues of the Dirac operator has 
already been discussed in Section II. Although it is most 
likely that this peak is a lattice artifact which disappears 
in the continuum limit, it seems useful to compare its 
shape with the one predicted by the model of dilute in- 
stanton - anti-instanton gas [7] . In this model instantons 
and anti-instantons interact at the scale of several lattice 
spacings rather than at a physical scale [7] , thus its suc- 
cess can serve as another argument that a strong excess 
of small eigenvalues at T > T c is a lattice artifact. 

In the model proposed in [7] each near-zero mode 
is associated either with an instanton or anti-instanton 
and the matrix elements of the Dirac operator be- 
tween the modes i and j are approximated by = 
ho exp (—dij/D), where ho is some energy scale, dij is the 
minimal distance in lattice units between the instanton 
associated with the mode i and the anti-instanton asso- 
ciated with the mode j and D is a typical scale of inter- 
actions between instantons and anti-instantons induced 
by fermionic fields. Matrix elements between the modes 
associated with two instantons or two anti-instantons are 
assumed to be zero. With ho being equal to the inverse 
lattice spacing, the optimal value of D for our data is 
D = 1.3 ± 0.1. Eigenvalue distributions obtained from 
this model at the optimal value of D are also plotted 
on Fig. 5. These distributions agree well with lattice 
data in the range A = 20 MeV . . . 100 MeV, while for 
A < 20 M eV the density of eigenvalues of lattice Dirac 
operator is significantly lower. 

Thus at high temperatures the distributions of eigen- 
values within the characteristic peak near A = are well 
described by the model used in [7], with the interaction 
scale between instantons being fixed in units of lattice 
spacing. This again suggests that the modes which con- 
stitute this peak should turn into exact zero modes in 
the continuum limit. This conclusion immediately rises 
the question about the model which could describe the 
spectrum of small eigenvalues of Dirac operator at not 
very high temperatures above T c in the continuum limit. 



IV. CONCLUSIONS 

In this work we have measured the spectrum of low- 
lying eigenvalues of overlap Dirac operator and the value 
of the chiral condensate in quenched SU (2) lattice gauge 
theory with tadpole-improved Symanzik action. Our re- 
sults suggest that depending on the temperature, the 
theory can be in the three regimes: a conventional con- 
finement phase at < T < T c w 320 MeV, for which 
the Polyakov line is zero and the chiral condensate is 
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FIG. 4: Distributions of low-lying eigenvalues of overlap Dirac operator compared with universal distributions for chiral 
orthogonal ensemble of random matrices. At the top: results for 10 3 x 4 lattice at /3 = 2.93, T = 264 MeV, at the bottom: 
results for 12 3 x 6 lattice at (3 — 3.23, T = 294 MeV, on the left: distributions of all eigenvalues, on the right: distributions of 
the lowest nonzero eigenvalues. 



nonzero, a phase at T > 1.5 T c » 480 MeV, for which 
the Polyakov line is nonzero and the chiral condensate 
is either very small or zero and which is therefore very 
similar to the conventional deconfinement phase, and a 
sort of "transition regime" at T c < T < 1.5 T c , where 
the Polyakov line is nonzero and the chiral condensate 
is still comparable with a typical hadron scale. The chi- 
ral condensate changes rapidly at the transitions between 
all these regimes, which suggests that even in quenched 
theory center and chiral symmetries arc closely related 
[17, 18], although in a somewhat more specific way than 
in the theory with dynamical fcrmions. 

The distributions of low-lying eigenvalues of the Dirac 
operator are universal in the confinement phase up to the 
critical temperature T c , and are described by the chiral 
orthogonal ensemble of random matrices, as it should be 
[2]. On the other hand, in the "transition regime" the 
distributions of small eigenvalues become non-universal. 
This could be expected, since the proof of the univer- 
sality of ps (EVA) and ps (T,VX min ) at small A is essen- 
tially based on the validity of the effective chiral the- 
ory [1, 2], although it is often claimed in the literature 



that the spectrum of small eigenvalues should be univer- 
sal whenever the eigenvalue density near zero is finite 
[2]. It is therefore reasonable to expect that at tem- 
peratures above T c , at which quarks are not confined 
in mesons anymore and the effective chiral theory be- 
comes invalid, these distributions become non-universal. 
In this temperature range the data are better described 
by the phenomenological model of dilute instanton - anti- 
instanton gas considered in [7] . It could be interesting to 
further investigate the applicability of such a model, since 
lattice simulations indicate that the density of topolog- 
ical charge is arranged into fractal-like structures which 
do not resemble instantons at all [19, 20, 21, 22]. Fur- 
thermore, a typical feature of the Dirac spectrum in the 
"transition regime" which is described by the model of 
7], namely, the characteristic peak near A = 0, seems to 
be a lattice artifact which disappears in the continuum 
limit. It is therefore an open question what kind of con- 
tinuum model may describe the properties of small Dirac 
eigenvalues above T c and what are the objects respon- 
sible for the persistence of near-zero eigenmodes of the 
Dirac operator. 
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FIG. 5: Eigenvalue distributions p (A) (on the left) and p(\ min ) (on the right) for 16 3 x 6 lattice at (3 = 3.275, T = 316 MeV 
compared with the phenomenological model of dilute instanton - anti-instanton gas used in [7]. 



Finally, at temperatures above approximately 1.5 T c 
the density of small eigenvalues rapidly drops down once 
more to a value which is very close to zero. Unfortu- 
nately, due to limited statistics it is very difficult to mea- 
sure this density and to establish whether it is exactly 
equal to zero or not. It is an intriguing open question, 
since for quenched theory chiral symmetry is not a sym- 
metry of the action at all and there are no theoretical 
reasons to expect that the chiral condensate turns to zero 
at some temperature. If it really does, this would mean 
that in pure gauge theory there are in fact two finite- 
temperature phase transitions, and it is not clear what 
might be the origin of the second transition. It is interest- 
ing to note here that some gauge theories constructed us- 
ing the gauge/gravity duality, such as the Sakai-Sugimoto 
model, indeed exhibit two finite-temperature phase tran- 
sitions, a deconfmcment one and the one associated with 



spontaneous breaking of chiral symmetry [23, 24]. 
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